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We investigate the general features of the renormalization-group flow at the Berezinskii-Kosterlitz- 
Thouless (BKT) transition, providing a thorough quantitative description of the asymptote critical 
behavior, including the multiplicative and subleading logarithmic corrections. For this purpose, we 
consider the RG flow of the sine-Gordon model around the renormalizable point which describes the 
BKT transition. We reduce the corresponding /3-functions to a universal canonical form, valid to all 
perturbative orders. Then, we determine the asymptotic solutions of the RG equations in various 
critical regimes: the infinite-volume critical behavior in the disordered phase, the finite-size scaling 
limit for homogeneous systems of finite size, and the trap-size scaling limit occurring in 2D bosonic 
particle systems trapped by an external space-dependent potential. 
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I. INTRODUCTION 

The Berezinskii-Kosterlitz-Thouless (BKT) theory 0- 
Al describes finite-temperature transitions in two- 
dimensional (2D) systems with a global U(l) symme- 
try, which belong to the so-called 2D XY universality 
class. BKT transitions are quite peculiar, since the low- 
temperature phase is not characterized by long-range 
order and the emergence of a nonvanishing order pa- 
rameter || Q, but rather by quasi-long range order 
(QLRO) with correlations decaying algebraically at large 
distance. For example, a 2D fluid of identical bosons 
cannot undergo Bose-Einstein condensation. Above T c 
these systems show a standard disordered phase with ex- 
ponentially decaying correlations. The BKT theory pre- 
dicts an exponential increase of the correlation length 
when approaching the transition point T c from above, as 
£ ~ exp(c/-y/r) with r = T/T c — 1. BKT transitions are 
generally expected in 2D systems of interacting bosonic 
atoms, such as those investigated in experiments with 
trapped atomic gases @-[Hj], in liquid helium films [ijj], 
in arrays of Josephson junctions [131 ]. etc. These experi- 
ments have provided evidence of the general features pre- 
dicted by the BKT theory. 

A standard representative model of the 2D XY univer- 
sality class is the classical 2D XY model. Its Hamiltonian 
is 



H 



XY 



i>i e u(i) 



(i) 



where the sum runs over the bonds of a square lat- 
tice. Its phase diagram shows a BKT transition be- 
tween a high-temperature disordered phase and a low- 
temperature QLRO phase, at [lj-yjj J/3 C = 1.1199(1). 
The asymptotic behaviors at the BKT transition, and in 
particular their logarithmic corrections, have been much 
investigated, see e.g. Refs. jla428j . by analytical and nu- 
merical studies. 

In this paper we investigate the general features of the 



renormalization-group (RG) flow at the BKT transition, 
providing a complete characterization of the asymptotic 
BKT behaviors, and in particular of the multiplicative 
and subleading logarithmic corrections. For this purpose, 
we exploit the mapping between the 2D XY or Coulomb- 
gas models and the sine-Gordon (SG) model, whose RG 
flow around the renormalizable point describes the BKT 
transition 34]. In order to investigate its RG flow, we 
first reduce the SG /3-functions to a canonical universal 
form. Explicitly, we show that we can define appropriate 
couplings u and v so that the associated /3-functions have 
the form j3 u = —uv and (3 V — — u 2 [l+vf(v 2 )] to all orders 
of the perturbative expansion in powers of u and v. The 
universal function f(v 2 ) cannot be determined by general 
arguments, but only by means of detailed calculations in 
the SG model (at present only f(Q) is known). Then, 
generalizing the RG study of Ref . [24j , we determine the 
asymptotic solutions of the RG equations in some differ- 
ent critical regimes: the infinite- volume critical behavior 
in the disordered phase, the finite-size scaling_[FSS) limit 
for homogeneous systems of finite size (35l l3a ] , and the 
trap-size scaling (TSS) limit [37j in 2D bosonic particle 
systems trapped by an external space-dependent poten- 
tial. The latter results are relevant for experimental in- 
vestigations of trapped quasi- 2D atomic gases, such as 
those reported in Refs. [7Hll| . 

The paper is organized as follows. In Sec.[TT]we summa- 
rize some of the RG ideas that we use to analyze the BKT 
RG flow. Sec. IIIII reports the derivation of the canonical 
form of the /3-functions and outlines the main features of 
the RG flow which they describe. In Sec. lIVI we derive the 
asymptotic critical behavior of some observables, such 
as the correlation function (■tjj(x)t/j(y)) of the XY model 
and its low-momentum components, when approaching 
the critical point from the high-temperature phase. In 
Sec. [V] we discuss the asymptotic behavior in the FSS 
limit, i.e. in the infinite- volume limit keeping the ratio 
between the system size and the correlation length fixed, 
and in particular at T c . Multiplicative logarithms appear 
in the two-point function and in the related susceptibil- 
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ity. Sec. lVIl is devoted to an analysis of the TSS behavior 
at the BKT transition. We show that, at T c , the criti- 
cal behavior with respect to the trap size must include 
new multiplicative logarithms. Finally, in Sec. IVIII we 
draw our conclusions and summarize the main results of 
the paper. The various appendices report some technical 
details of the derivations of the results. 



Therefore the flow line lies in the intersection of all the 
Hi hypcrsurfaces. 

We can also take into account the size L of the system. 
With the usual hypotheses [39M4l| of the FSS theory, this 
is achieved by adding a term L/b in the scaling Ansatz: 

F(u w ,u 20 ,...,L)=b- d F[u 1 (b),u 2 (b),...,L/b}. (9) 



II. A SHORT SUMMARY OF THE 
RENORMALIZATION-GROUP IDEAS 

Before discussing the RG flow at the BKT transition, 
we would like to recall a number of ideas concerning the 
RG flow, which we then apply to the study of the BKT 
critical behavior. We consider a generic critical system 
depending on a set of scaling fields u\o, u 2 q, and so on. 
The RG flow with respect to a length rescaling b is defined 
by 1H 



(2) 



with boundary condition Ui(b = 1) = u i0 . Then, the 
scaling part of the free-energy density satisfies 



J r (u w ,u 20 ,---) = b F[ui(b),u 2 {b), 



(3) 



where d is the space dimension. Analogously, the corre- 
lation length satisfies 

£("10, "20, • • •) = b£[ui(b), u 2 (b), ...}. (4) 

An operator that renormalizes multiplicatively scales as 

O(uio,uz , ...) = (5) 
= b d °Z [ Ul (b), u 2 {b), . . ]0[ux{b), u 2 (b), . . .], 

where do is its power-counting dimension and Zq satis- 
fies the RG equation 



d In Zo 
} ~db~ 



7o0i,w 2 , • ■ •): z o(b = !) = !, (6) 



where 70 is the anomalous dimension associated with O. 

The RG flow can be characterized in terms of char- 
acteristic surfaces. For a given set of initial conditions 
we consider the functions fi(u\) for i > 2 (the choice 
of the first scaling field to parametrize the flow is arbi- 
trary; any other choice would work equally well), which 
are solutions of the equations 



dfi _ Pi[ui, f 2 {ui), f 3 (ui), 
dui f3i[ui, f 2 (ui), f 3 {ui), 



(7) 



with initial conditions /i(uio) = um- The hypersurface 
Hi of equation Ui = fi{u\) is invariant under the RG 
flow, since 



^[«i(i)-A(«i(6))]=A-|A 



0. 



(8) 



Also correlation functions that depend on coordinates x, 
y, etc. can be analyzed. In this case the RG mapping is 
simply x — > x/b, y — > y/b, etc. 

In the following we shall use the notation 



I = ln6, 



- = b- 
dl db' 



(10) 



III. 



CANONICAL FORM OF THE BKT BETA 

FUNCTIONS AND 
RENORMALIZATION-GROUP FLOW 



We want to study the general features of the RG flow 
at the BKT transition of 2D systems with U(l) symme- 
try. For this purpose, we consider the SG field-theoretical 
model, see, e.g., Ref. [42|, defined by the Lagrangian 



C 



SG 



a 2 f3- 



[1 - cos(m 



(11) 



where a and (3 are dimensionless coupling constants, and 
a is an ultraviolet length scale. The RG flow around the 
fixed point /?* = \/8~7r, a* = describes the BKT transi- 
tion [34| . It can be investigated by a renormalizable two- 
parameter perturbative expansion in powers of 6, defined 
by (3 2 = 8tt(1 + S), and a. Their j3 functions have been 
computed to two-loop order [§3, 51] , obtaining Q 



P a = 



-2aS-—a 3 , 
64 ' 

1 2 1 2r 

32 a + 16 aS - 



(12) 
(13) 



Under an appropriate analytic redefinition of the cou- 
plings a and 6, the above two-loop /3-functions can be 
rewritten as 



-UV. 



1 2 

-u u v. 

2 



(14) 



The coefficient —3/2 is universal in the following sense: 
there is no redefinition of the couplings v! = U(u,v) and 
V(u,v) such that /3 U < = —u'v' and = — u' 2 — 
+ . . ., with c ^ -3/2. 
This reduction to a universal form can be extended to 
all orders, leading to the most general canonical form of 
the /3-functions which is compatible with the invariance 
of the model under a — ¥ —a (it corresponds to a shift 



/2 
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of ir/p in the field <f). We prove that, by an analytic 
redefinition of the couplings 



n+m>2 



^2 a S,n m U n V m , 



(15) 
(16) 



n+m>2 



the /3-functions of the SG model can be rewritten in the 
general form 



Pv{u,v) 



-uv, 

,2r 



(17) 

(18) 



-u*[l + vf(v% 

where f(x) has an expansion of the form 

f(x) = b Q + b 1 x + b 2 x 2 + ... (19) 

This representation of the /^-functions is universal, in the 
sense that, by redefining the couplings, it is not possible 
to obtain f3 functions of the same form, i.e., (3 u i = —u'v' 
and (3 v i = — u' 2 [l+v'g(v' 2 )], with g(x) ^ f(x). The proof 
is outlined in App.[X](this result was already conjectured 
in Ref. [45;] without proof). The coefficient bo can be read 
off from Eq. ([Til), 



3/2, 



(20) 



while the higher-order terms, bi, i > 1, are unknown. 
Note that the evaluation of the next universal coeffi- 
cient bi requires a perturbative calculation of the SG 
/3-functions to four loops. 

The analysis of the RG flow can be performed following 
the method outlined in Sec. HH see also Refs. [H, |34| . 
First, we define the RG invariant function 



where 



F(v) 



Q(u,v) = u 2 -F(v), 



wdw 



o l + wf(w 2 ) 



(21) 



(22) 



which satisfies 



dQ dQ ^dQ 



0, 



(23) 



where I is the flow parameter. The RG flow follows the 
lines Q = constant. It is thus natural to parametrize the 
RG flow in terms of Q and v(l). Since 

^ = &(«,«) =-[Q + F(v)][l + vf(v 2 )}, (24) 



we obtain 



/ = - 



v(l) 



dw 



[Q + F(w)][l+wf(w 2 )Y 



(25) 



where v(l = 0) = vo- 
lt is important to stress that Eqs. ([TT|) and ([T%f are 
the two /3-functions associated with the marginal oper- 
ators that characterize the BKT transition. For a full 
understanding of the scaling corrections one should also 
consider the contributions of the subleading operators, 
which are suppressed by powers of the critical length 
scale. The most relevant subleading operator at the BKT 
transition is expected to have RG dimension —2, as in 
the Gaussian spin- wave theory, see, e.g., Ref. 1141. In 
the standard RG language this corresponds to a scaling- 
correction exponent u) = 2. Thus, subleading operators 
induce corrections of order £ -2 in the high-temperature 
infinite-volume limit and of order L~ 2 in the FSS limit 
(apart from corrections arising from the boundary con- 
ditions 5SE3i which are expected to be 0(L _1 ); they 
are absent in the case of boundary conditions preserv- 
ing translation invariance, such as periodic boundary 
conditions). Additional multiplicative logarithms may 
also appear (hence, corrections might scale as £ -2 ln p £, 
L~ 2 \n q L), because of the possible resonance between the 
subleading and the marginal operators [381 ) , the difference 
between their RG dimensions being an integer number. 
In the following we shall not consider these scaling cor- 
rections, since our focus will be mainly on the logarithmic 
corrections to the leading behavior that can be obtained 
by considering only the two marginal couplings. 



IV. INFINITE- VOLUME RESULTS AT THE 
BKT TRANSITION 

Let us now apply these results to the XY model. In 
Fig. [1] we show a sketch of the RG flow. Repeating the 
discussion of Refs. [!, IHHl] the XY model can be mapped 
onto a line in the (u, v) plane with v > 0. The BKT 
transition is the intersection of this line with the line 
Q = and the high-temperature phase corresponds to 
Q > (region C with v > in Fig. [I}. Thus, Q plays 
the role of thermal nonlinear scaling field, i.e. 



Q = qiT + qiT 



(26) 



where r = T/T c — 1, and qi are nonuniversal coefficients. 
In the following we shall use Q instead of r. 

Let us first consider the infinite-volume correlation 
length £oo(t) = £,oc(Q,v), which we may define by using 
the second-moment of the two-point correlation function 



G(x,y) = (iKxMy)), 



(27) 



or its large-distance exponential decay. For all v > the 
correlation length is singular as Q — > 0. On the other 
hand, for v < 0, the correlation length is analytic on the 
line Q = (see Ref. [24j]). Hence, in order to obtain the 
singular behavior, we use the RG equations to flow from 
the starting point vq > to the negative point v = — 1. 
We determine Iq by requiring 



v(l ) 



(28) 
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FIG. 1: Sketch of the RG flow at the BKT transition. The 
dashed curve in the region C (v > 0) shows the approach to 
criticality (Q = line) of the 2D XY model from the HT 
phase. The correlation length is singular along the Q = line 
in the region v > 0, while it is analytical along the Q = for 
v < 0. 



so that 



rfi/> 



J(Q,«o). (29) 



Then, Eq. (g]) gives 

£oo(Q,«o) =e'£ oo (Q,u(0) = e / «^ 0o (Q ) -l). (30) 

Since £oo(Q, — 1) is analytic in Q, the singular part is 
given by the exponential term. The behavior of £oo(t) 
for r — > is obtained by computing the asymptotic ex- 
pansion of I(Q,vq) for Q — > 0, which can be written in 
the form 

I(Q, v ) = -^=Y, 7 "Q" + E Y nMQ n , (31) 

see App. [B] for its derivation. The nonanalytic terms in 
this expansion depend only of the coefficients bk which 
appear in Eq. (|Tg)) . The first two coefficients are 



h = 7T, 
h = 



12 



3tt 
16' 



(32) 



Correspondingly, we obtain 

Coo(r) =Xexp( 7 r/v / Q)[l + ^V / Q + 0(Q)]. (33) 

Expanding Q in powers of r we obtain the celebrated 
BKT expression for the correlation length 0, 9. Of 
course, this is also consistent with the exact result for 
the correlation length of the F-model, computed using 
transfer- matrix techniques, see e.g. [Til [29j. 

These results allow us to derive the critical behavior of 
of generic RG invariant dimensionless quantities, which 
we generically denote by R. Examples of such quan- 
tities are ratios of length scales and the so-called zero- 
momentum 2n-point renormalized couplings (related to 



the small-magnetixation expansion of the scaling free en- 
ergy), see e.g. Refs. [20I [30h331 | . Indeed, R satisfies the 
scaling relation 



R(Q,v ) = R[Q,v(l)}, 



(34) 



for any I. This implies that R(Q,vq) is independent of 
v(l), hence an analytic function of Q and therefore of r. 
It follows 



R(t) = R* 



cr 



+ O(ln- 3 £ OC ), (35) 



where the costant cr is expected to be universal. 

Let us now consider a generic operator that renormal- 
izes multiplicatively. Writing the corresponding anoma- 
lous dimension 70 (u, v) in terms of Q and v, from Eq. © 
we obtain 



Zo(Q,vo) 



, 1 , lo(Q,w) 
exp / aw — — — r 
' ' Pv{Q,w) 



Z (Q,-l). (36) 



Taking into account the symmetry properties of the SG 
model, we write the expansion of 70 as 

7o = .9oo + goiv + go2V 2 + g2QU 2 + 0{v 3 ,u 2 v 2 ,u 4 ) = 
= 9oo + 9oiv + (go2 + 92o)v 2 + g 2 oQ + 0(v 3 ,Qv 2 , Q 2 ). 

(37) 

It is important to discuss the universality of this expan- 
sion. There is a residual transformation of the couplings 
that leaves the /3-functions ([TT]) and (|T5)) invariant: 



U 
v' 



■ u ■ 
v - 



- Auv + 
Au 2 + . 



(38) 
(39) 



for any A. The RG invariant function Q is invariant 
under the transformation and so are the coefficients goo, 
goi, and 1702, hence they are universal. Instead, 320 can 
be changed at will, hence it is model dependent. 
We can now rewrite Eq. ([36)) as 



Zo(Q,vo 

„'oSoo 



(40) 



, 1 , jo(Q,w) 

exp / aw — — r 

1 ; Pv(Q,w) 



goo 



z (Q,-i). 



Collecting everything together we obtain 

Coo (Q, v) d °+9°° = ^ P J Vo dw M Q,w) > 

(41) 

where C{Q) is an analytic function of Q. For Q — > we 
obtain an expansion of the form 

0(Q,v) 



Coo(Q,w) do+9o ° 



aoC(Q)[l + a iy /Q + a2Q + ...], (42) 



where ao and 02 depend on nonuniversal details (the 
starting point vq, for instance), while a\ is universal, since 
it only depends on the universal coefficient go2'- 



a 1 



-7TS02- 



(43) 
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The above result can be specialized to the susceptibil- 
ity, defined as the space integral of the two-point function 
([27f . Perturbation theory for the scaling dimension of the 
spin correlation function gives \['> ij 



7= -- + -S- -<T + /iia 2 
AAA 



(44) 



where /ii is an unknown coefficient. If we perform the 
redefinitions (a, (5) — > (u, v) considered before, we can 
rewrite 7 as 



1 1 1 2 2 

1 = — A + %v - -v- + g 20 u . 



(45) 



with arbitrary (720- The previous results show that, in 
the infinite-volume limit, the susceptibility satisfies the 
scaling law 



xC 7/4 = A(l+ciyQ + c 2 Q+...). 



(46) 



The coefficient c\ can be computed exactly. Since (702 = 
— 1/16, we obtain 



Cl= 16- 



Using Eq. (J33j) we can write 



Q 



and obtain 



161n(W*) 



O(ln- 2 $oo) 



0(1/ In 2 Co 



(47) 



(48) 



(49) 



As already noted in Ref. [24]], RG predicts the absence 
of a leading logarithmic singular term. Note also that, 
at variance with what happens with RG invariant quan- 
tities, cf. Eq. (1351) . corrections decay as l/ln^x, and not 
as 1/ In 2 Coo- Moreover, the coefficient of the leading log- 
arithmic scaling correction is universal. 

We should note that in the literature, see, e.g., 
EH, HI)], it was often assumed that the correct scaling 

behavior is xCoo 7 ^ 4 ~ (In Coo) 1 / 8 - Such a behavior was 
derived as follows. Since the two-point function at the 
critical point scales as [US El] G(r ) ~ (hir) 1 / 8 /r 1 / 4 , it 
was argued that 



X 



d 2 rG{r) ~ ^(Infoo) 1 ' 



(50) 



However, in the presence of logarithmic corrections it is 
not clear if one should cutoff the integral at Coo , or rather 
at Coo (In Coo) 9 - If we use the latter cutoff we can freely 
change the leading logarithmic behavior by changing q. 
In particular, if q = —1/14, we obtain Eq. (|4"5|) . The most 
recent MC simulations of the XY model in the thermo- 
dynamic limit confirm the absence of the leading singular 
term [251 ] . The analyses of high-temperature expansions 
have led to apparently contradictory results: the results 



Ref. [27l| ar e consistent with Eq. (|49|) . while the analysis 
of Ref. [28l ] apparently favors the presence of a leading 
logarithmic correction. On the other hand, we should 
note that the resummation of high-temperature expan- 
sions at finite order can hardly reconstruct an asymp- 
totic behavior at the critical point with logarithmic cor- 
rections. 



V. FINITE-SIZE SCALING 

A. Finite-size scaling in the high-temperature 
phase 

In order to study the FSS regime, we add L/b = Le~ l 
in the scaling Ansatz. If Q ^ 0, i.e. we are not at the 
critical point, we can study the FSS regime as we did 
in the previous section. If we choose I = lo such that 
v(Iq) = —1, we can write 



Le~ l ° =Le- I ( Q > v °) =&o(Q,-l) 



L 



(51) 



Coo (Q, Wo)' 

where Coo(Q, vo) is the infinite- volume correlation length, 
and Coo {Q, — 1) is an analytic function of Q, which is finite 
for Q — !• 0. The finite-size correlation length must satisfy 
the equation 

C(Q,«o,i) = e i »C(Q,-l,ie- i «). (52) 
For Q — > 0, introducing the FSS variable 



we obtain 

£,(Q,v ,L) 
Coo (Q, v a ) 



z — Coo/ L, 



£(Q,-l,Le- 



Cc 

A(z) 



+ QB(z) 



(53) 



(54) 



0(Q 2 



Hence, if we use z as basic FSS variable, scaling correc- 
tions decay as 1/ In 2 Coo ■ It is clear that the same argu- 
ments apply to any observable, as long as we divide it by 
its infinite- volume limit. Hence 



O{Q,v ,L) <D(Q,-l,e- l °L) 



e»oo(Q,-i) 

A (z) + QB (z) 



(55) 



0(Q 2 ). 



All these relations hold as long as z is finite. The infinite- 
volume limit is not uniform in z and indeed the scaling 
functions Ao(z), Bq(z), etc, are singular for z — > 00. 
The limiting behavior for z — > 00, i.e. the asymptotic 
behavior at the critical point will be discussed in the 
next section. 



B. Finite-size behavior at the critical point 

The finite-size behavior at T c is not simply obtained 
by extending the results of Sec. IV Al to T c . We consider 
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the correlation length in a finite box of size L, £(Q, v, L). 
For finite L, we can take the limit Q — > and obtain 
£(0, v, L), which is singular as L — > oo. At T = T c it is 
convenient to fix I = In L, which gives 

£(0, v, L) = eJ^O, v(l), e- l L] = L£[0, v(L), 1], (56) 

where v(L) — v(l — InL). The function £[0,v(L),l] is 
analytic for any L, since the size is equal to 1, hence we 
have 

£(0, v, L) = L[s + S!V(L) + s 2 v(L) 2 + ...], (57) 

where the coefficients Si are universal and v(L) is defined 
by 



InL 



I'll 



dw 



V(L) F(w)[l + wf(w*)Y 



(58) 



For L — > oo, the effective coupling v(L) vanishes as 
1/lni. In this limit we obtain 



lni 



b 



v(L) 



lnv(L)+K 



(59) 



v(L) 



dw 



1 L + ^L 

F(w)[l + wf(w 2 )} w 2 3w 



where 



K = / dw 

Jo {F(w)[l+wf(w 2 )} w 2 ' 3w 



A + (60) 



In Eq. (|59l) all terms are universal, except for the constant 
A, which encodes all microscopic details. If we define 



H = ln(Le~ K ), 
for v{L) — > we obtain the expansion 



(61) 



1 h ^h 2 

M = ^ + j + -^v(L) + 0(v 2 ). (62) 



The general solution is 



v(L) 



In fi In ji 



- + CTi^ r + 0'2- , 



(T 3 - 



ln/z 



+ 



C4^t + 0(m 4 In 3 /i) 



(63) 



where all coefficients of the expansion are universal, and 
<Ji up to i = 4 can be computed in terms of bo only. The 
terms of order In" -1 can be resummed, obtaining 



v(L) 



1 



[i + \ In 



+ 0(^ 3 ln/z). 



(64) 



Indeed, we can rewrite Eq. (|fj2"j) as 



— ~/i ^ In In /i = /i+ - In /i. (65) 



The above results allow us to derive the asymptotic finite- 
size behavior at T c of generic RG invariant dimension- 
less quantities R, such as ratios £/L for any definition of 
length scale, Binder cumulants and the helicity modulus 
T. They are expected to behave as 



R{L) = R* + C R v(L) + 0(v 2 ) 



(66) 



where the R* and Cr are universal, although they may 
depend on the shape of the finite volume and the bound- 
ary conditions. Then, using Eq. (|64[) , we obtain 



i?« R* 



C R 



/* - 



\ ln^ 



0(/T 



(67) 



This result improves earlier asymptotic expansions, see, 
e.g., Refs. [15L l48l450| . The asymptotic values R* and Cr 
can be c omp uted within the spin-wave theory, as shown 
in Refs. [l5l |48| . For example, in the case of the helic- 
ity modulus in a square lattice with periodic boundary 
conditions, T* = 0.636508... and C T = 0.318899... [l]f. 
An analogous result applies to the exponential correla- 
tion length £ e in a strip L x oo, i.e. for i?^ c = £ e /L, with 
Rt = 4/tt and Q e = 2/tt 0. 

It is worth noting that the behavior of RG invariant 
quantities at T c has different features with respect to the 
infinite- volume case. Here corrections decay as 1/lnL, 
while in the infinite- volume case, see Eq. (|35p , they decay 
as 1/ (ln^oo) 2 , i.e. with the square of the logarithm of the 
relevant length scale. 

The finite-size behavior of observables O with anoma- 
lous RG dimension can be obtained in a similar fashion. 
We need to compute the behavior of the integral 



dw- 



7o(0,w) 



/M0,™) Ml) F(w)[l+wf(w 2 )} 
7o(0,uO - .goo 

dw- 



= goo InL - ffoi Int'(^) + K' + 0(v), 



(68) 



where 70 (Q, w) is the anomalous dimension as a function 
of Q and v, and A'' is a nonuniversal constant. Hence 

Z[0,v(L)} = L 900 v(L)- 901 e K ' [1 + O(v)]Z(0,v ). (69) 
We end up with 

The above results imply that the two-point function at 
T c behaves as 

/ 1 \ 1/8 
G(x,y)«L- 1 /M M+ _ln / iJ Q( X /L,y/L). (71) 

In particular, its space integral, i.e. the susceptibility, 
scales at the critical point as 

XL- 7 '* = A L + I htfj) V8 [1 + OO*" 1 )] , (72) 
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where we used goi = 1/8. A numerical analysis of the 
2D XY model providing evidence of the lea ding multi- 
plicative logarithm is reported in Refs. [l5|, Il6j . Note 
that a naive integration of the infinite- volume two-point 
function G{r) up to r ~ L would give the same re- 
sult, x ~ i 7/4 (lni) 1/8 . The critical-point behavior ([72]) 
should be contrasted with Eq. (|49[) : in infinite- volume 

X ~ Coo 4 without additional leading logarithms. 

VI. TRAP-SIZE SCALING 

Statistical systems are generally inhomogeneous in na- 
ture, while homogeneous systems are often an ideal limit 
of experimental conditions. Thus, in the study of critical 
phenomena, an important issue is how critical behaviors 
develop in inhomogeneous systems. Particularly interest- 
ing physical systems are interacting particles constrained 
within a limited region of space by an external force. This 
is a common feature of recent experiments with diluted 
atomic vapors [H| and cold atoms in optical lattices [52j , 
which have provided a great opportunity to investigate 
the interplay between quantum and statistical behaviors 
in particle systems. 

Experimental evidences of BKT transitions in trapped 
quasi- 2D atomic gases have been reported in Refs. 0- 
Hl| . The inhomogeneity due to the trapping potential 
drastically changes, even qualitatively, the general fea- 
tures of the critical behavior. For example, the corre- 
lation functions of the critical modes do not develop a 
diverging length scale in a trap. Nevertheless, when the 
trap size becomes large the system develops a critical 
scaling behavior, which can be described in the frame- 
work of the TSS theory [§3, [H| . TSS has some analogies 
with the standard FSS for homogeneous systems with two 
main differences: the inhomogeneity due to the space- 
dependence of the external field, and a nontrivial depen- 
dence of the correlation length on the trap size at the 
critical point. 

The above considerations apply to general quasi-2D 
systems of interacting bosonic particles trapped by an 
external harmonic potential. In particular, we mention 
systems of bosonic cold atoms in quasi- 2D optical lat- 
tices [Hf, which can be effectively described [54| by the 
Bose-Hubbard (BH) model (H 

#bh = -^(blbj + bth) (73) 

{ij) 

+-^^2ni(m - l) - n^m, 

i i 

where hi is the bosonic operator, rii = b\bi is the particle 
density operator, and the sums run over the bonds {ij) 
and the sites i of a square lattice. The phase diagram of 
the 2D BH model presents finite-temperature BKT tran- 
sition lines, connecting T — quantum transitions, such 
as those from the vacuum state to the superfluid phase, 



and from the superfluid phase to Mott phases [H|. Ex- 
periments with cold atoms (5ll l52j are usually performed 
in the presence of a trapping potential, which can be 
taken into account by adding a corresponding term in 
the Hamiltonian: 

Ht^ = H^ + Y j y{n)n l , (74) 

i 

V(r) = w p r p , (75) 

where r is the distance from the center of the trap, and p 
is a positive even exponent. A natural definition of trap 
size is provided by 

l t = J^ p /w. (76) 

Far from the origin the potential V(r) diverges, therefore 
(rij) vanishes and the particles are trapped. The trapping 
potential is effectively harmonic in most experiments, i.e. 
V = 2. 

The trapped 2D BH model has been numerically in- 
vestigated in Ref. [So ], by quantum Monte Carlo simu- 
lations, showing that the BKT critical behavior is sig- 
nificantly modified by the presence of the trap. Anal- 
ogously, an accurate experimental determination of the 
critical parameters, such as the critical temperature, crit- 
ical exponents, etc., in trapped particle systems requires 
a quantitative analysis of the trap effects. In the follow- 
ing we investigate this issue at the BKT transition and 
derive the asymptotic TSS behavior from the BKT RG 
flow. 



A. General features of trap-size scaling 

Let us first describe the general TSS approach to stan- 
dard continuous transitions [33| , characterized by two rel- 
evant parameters r and h (usually r ~ T/T c — 1 and h is 
the external field coupled to the order parameter) , whose 
RG dimensions are y T — 1/v and yh = (d + 2 — rj)/2. 
The presence of an external space-dependent field V(r) = 
(wr) p significantly affects the critical modes, introducing 
another length scale, the trap size k ~ 1/w. Within the 
TSS framework [§3, [Hj], the scaling of the singular part 
of the free-energy density around the center of the trap 
is generally written as 

F(x,T,h) = r t ed T(rl-\rf t ^ ,hl 6 t y % (77) 

where r is the distance from the center of the trap, and 
9 is the trap exponent. TSS implies that at the criti- 
cal point (r = 0) the correlation length £ of the critical 
modes is finite, but increases as £ ~ l\ with increasing the 
trap size It- TSS equations can be derived for the corre- 
lation functions of the critical modes. For example, the 
correlation function of the fundamental field i[>(x) (the 
quantum field b in the BH model) is expected to behave 
as 

G(x,y) = $(x)V>(y)) c = i; ev G(^ 9 ,ylt e ,rf t /u ), 

(78) 
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where Q is a scaling function. 

The trap exponent 9 generally depends on the univer- 
sality class of the transition, on its space dependence (in 
experiments the external potential is usually harmonic), 
and on the way it couples to the particles. Its value can 
be inferred by a RG analysis of the perturbation Py rep- 
resenting the external trapping potential coupled to the 
particle density. The universality class of the superfluid 
transition can be represented by a $ 4 theory for a com- 
plex field ip associated with the order parameter, see, e.g., 
Ref. El, 



taking into account the trapping potential shows that 
the asymptotic trap-size dependence at the BKT tran- 
sition presents multiplicative logarithms which turn out 
to depend on the power law of the trapping potential. 
Therefore, at a BKT transition the TSS relations ((77)) 
and (|78[) must be revised, including multiplicative loga- 
rithms, which differ from those observed in the FSS case. 



B. Renormalization-group analysis of trap-size 
scaling 



H, 



<i> 1 



d d x [|c^(x)| 2 + r#(x)| 2 + W |V(x)| 4 ] . (79) 



Since the particle density corresponds to the energy op- 
erator \ip\ , we write the perturbation Py as 



P v = / d d xV(x)\i/>(x)\ s 



(80) 



The exponent 6 is related to the RG dimension y w of the 
coupling w of the external field V = (wr) p by 9 = l/y w . 
Then, a standard RG argument gives 



py w -p + Vn 



(81) 



where y n = d — \jv is the RG dimension of the den- 
sity/energy operator l^l 2 . We eventually obtain 



1 

Vw 



pv 



1 



pv 



(82) 



We may derive the value of 9 at the BKT transition by 
formally setting v = oo in Eq. (|82j) . somehow correspond- 
ing to the BKT exponential behavior of the correlation 
length £ - exp^r" 1 / 2 ), where r = T/T c - 1 -> Q+. This 
would give 9 = 1 for any power p of the potential. This 
result is also obtained by extending to the BKT transition 
point the result 9 = 1 for the TSS in the whole QLRO 
phase [57[, which can be inferred by a RG analysis of 
the trap perturbation along the low-temperature line of 
Gaussian fixed points where spin-wave theory applies. 

The trap-size dependence predicted by TSS has been 
verified at various phase transitions, for example at the 
Ising transition of lattice gas models [13, [HI , at the 3D 
superfluid transition of bosonic particle systems such as 
those described by the 3D BH model (59|, and at the 
quantum T = Ising and Mott transitions [H, I60l - l62j | . 
Multiplicative logarithms are generally expected at the 
upper dimension of the given universality class. We shall 
show that they also appear at the BKT transition, which 
should not be surprising because they are already present 
in the scaling behavior of homogeneous systems, as dis- 
cussed in the previous sections. 

Note that the RG dimension associated with the size 
L is also 1 (in length units). This might suggest that 
TSS is analogous to FSS, i.e. characterized by the same 
power laws and similar multiplicative logarithms. How- 
ever, as we shall see below, the analysis of the RG flow 



To investigate the TSS regime, we extend the RG anal- 
ysis presented above. It is quite obvious that the presence 
of the trap does not change the short-distance behavior 
of the model, hence no change should be made on the 
scaling behavior of the couplings. Let us now consider 
the flow of the coupling w entering the potential (|74[) . 
which, in full generality, can have the form 



dw 

~dJ 



fi w (u,v,w). 



(83) 



If we start from w = (no trap), we should always have 
w(l) = 0, hence the /3-function should have the form 

p w {u, v, w) = wH (u, v, w). (84) 

Assuming y w = 9 = 1, we have 

H(0,0,0) = y w = l. (85) 

In App.Owe show that we can define a nonlinear scaling 
field z(u,v,w) so that the /3-function (|84|) becomes 



/3 z (u,v,z) = zT(Q,v), 



(86) 



where, as before, we have replaced u with the RG invari- 
ant quantity Q and T(0,0) = 1. The RG flow of z is 
particularly simple: 



z(l) = z exp{ / T[Q lV {l')}dl 




(87) 



r {l) , T(Q,w)-l 1 

Z ° e ' eXP '-l V0 dw lQ + FM]H + W f^)]j> 

where zq ~ 1/h is the starting point of the flow. Below, 
we shall consider two cases: first, we consider the high- 
temperature phase T > T c , then TSS at the critical point. 
In both cases, we assume that the infinite-volume limit 
has been attained, i.e. that L 3> l t . 



1. TSS in the high-temperature phase 

We start from the general scaling relation for the two- 
point correlation function 

G(x, y; Q, v , z ) = Z[Q, v(l)]G[x.e~ l ,ye~ l ,Q, v(l), z(l)]. 

(88) 
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Note that the renormalization function Z(Q, v) does not 
depend on the scaling field z, since the renormalization 
constant is only determined by the short-distance behav- 
ior of the operators — the fundamental field in this case 
— defining the correlation function. As in Sec. lIVl we fix 
I = lo by requiring v(lo) = — 1. Then, by using Eq. (|30|) 
we can write 

e'° = 5°°^' ^ w a^(Q, vo)[l + O(Q)], (89) 

where a = l/£oo(<3 = 0,-1) is a constant. To obtain the 
trap corrections, we must evaluate z(Iq). Eq. (|57|) allows 
us to write 



z('o) = -z a^ 00 ((3,i;o) 



x exp 



dw 



T(Q,w)-l 



[Q + F(w)][l + wf{w 2 )] 



The integral is finite for Q — > 0, with corrections of order 
y/Q, see App. [Bj Hence 

^ao)^^oo(g,« )[i + o(Q 1 / 2 )]. (91) 

Since It ~ l/^o, we obtain the general scaling form 

G(x,y;r) -g(xC 1 ,yC 1 ,^C 1 )+0( v / Q). (92) 

2. T55' at the critical point 

Let us now consider TSS at criticality {Q = 0). In this 
case it is convenient to fix Z = 1$ by setting z(lo) = 1. 
For Z — > co, v(l) —> 0, hence we can rewrite Eq. (|87| as 

l = ^Toe- io ua o )-* 1 [l + 0(^(/o))], 

where ii is defined by the expansion 

T(0,w) = 1 + tiW + Oiw 2 ), 

and Tq is a nonuniversal constant given by 



(93) 
(94) 



Tq = v li exp 



dw 



T(0,w) - 1 _ h 
F(w)[l + wf(w 2 )} ~ w 



(95) 

The flow of v(l) follows from Eq. (|25|) . Using the results 
of Sec. IV Bl and in particular Eq. (f5T?|). we have v(l) = 
l^ 1 + 0(l~ 2 In I), so we obtain 

l = z T e l l tl [l + O(r 1 lnl)]. 

Inverting this equation we obtain 

i IhizoTol-* 1 



(96) 



zqTq 



(97) 



To obtain the scaling behavior, we need to determine the 
large-Z behavior of Z[0, v(l)]. Using the results of Sec. fVBl 
and in particular Eq. (|69[) . we obtain 



Z[0,v(l)} 



Jgo 



v(l)- gm e K ' Z(0,v ) 



e z(0,u ). (98) 



Substituting this result into Eq. (|88[) . and choosing the 
length scale so that zqTq w + 0(/ t -2 ), we obtain the 
TSS of the two-point function at T c . We write it as (using 
ffoo = -1/4 and 501 = 1/8) 

G(x,y) = / t - 1/4 (lnZ,) 1/8+K/4 x (99) 
x g[ X (\nl t ) K /l t ,y(lnl t r/l t }, 

where k = t\. 

We may also define the trap susceptibility Xt, 



XL 



(90) and the trap correlation length £ t , 



Eq. ([99]) implies the asymptotic behaviors 



^ ~ lt(\nl t )' K . 



(100) 



(101) 



(102) 
(103) 



We do not know the value of the coefficient t\ appear- 
ing in the expansion (|94j) . which provides the exponent k 
in the asymptotic formulas. We generally expect that it 
depends on the power p of the trapping potential. Note 
that k must vanish in the limit p — > 00. Indeed, in this 
limit the trapped system is equivalent to a homogeneous 
system confined in a circle of radius It — 1/w with open 
boundary conditions. Therefore, TSS is equivalent to 
standard FSS with L ~ l t , hence k = 0. We anticipate 
that the numerical analysis that we present below pro- 
vides a strong evidence that k depends onp; in particular, 
it suggests k = 2/p. 



C. Numerical results in the presence of an external 
space-dependent field 

The main features of TSS are expected to be univer- 
sal, hence the RG results should apply to generic 2D 
systems characterized by a U(l) symmetry in the pres- 
ence of an external space-dependent field coupled to the 
energy density. For a numerical check of the RG pre- 
dictions, we consider the classical 2D XY model in the 
presence of an external space-dependent field coupled to 
the energy density. The Hamiltonian is given by 

H v = - jJ^Refe^i, (104) 

m 

Uij = [1 + V{ nj )]-\ V(r) = W P r P, (105) 



r.ij is the distance 



where p is an even positive integer 
from the origin of the midpoint of the lattice link con- 
necting the nearest-neighbor sites i and j. We set J = 1. 
The inhomogeneity arising from the space dependence of 
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FIG. 2: FSS plot of the two-point function G(0, x) at T c for 
x = (r, 0), < r < k — (L — l)/2, for a homogeneous XY 
model on a L 2 lattice with open boundary conditions (OBC). 



Uij is analogous to that arising from a trapping potential 
in particle systems. Thus, It = 1/w may be considered 
as the analog of the trap size (f76|) . For p — >• oo, V(r) =0 
for r < l t and V(r) = oo for r > l t . Hence, the system is 
equivalent to a homogeneous system confined in a circle 
of radius It with open boundary conditions. Therefore, in 
this limit TSS must reproduce standard FSS. A study of 
the tra p eff ects in the low-temperature phase is reported 
m Ref. [13. Here we focus on the trap-size dependence at 
the BKT critical temperature T c = 0.89294(8) [I1GJ] of 
the homogeneous XY model (JlJ , which is also the model 
(fT04f with Uij = 1. 

We present results of Monte Carlo (MC) simulations of 
model (|104p . We use a mixture of Metropolis and overre- 
laxation updates of the spin variables [50j • We consider- 
square lattices with L 2 , odd L, sites and open boundary 
conditions. Lattice points have coordinates (x, y) with 
—{L - l)/2 < x,y < (L - l)/2, so that the origin (0,0) 
is at the center of the lattice. The external potential is 
given by V(x) = (r/l t ) p , where r = |x| and l t is the trap 
size. The lattice size L is taken sufficiently large to avoid 
finite-size effects. We check that they are negligible com- 
pared with the statistical errors by comparing results at 
fixed trap size It for different lattice sizes L. 

We want to check the TSS equation for the corre- 
lation function at T c . For this purpose we report results 
at T c for the correlation function G(0,x) = (ip(0)ip(x)}, 
which is expected to scale as 

G(0,x) = l^^ilnlt) 1 ^ 4 g p [r(lnl t r/l t ], (106) 

where r = |x|, and the exponent k is expected to depend 
on the power of the external potential. 

We first check the FSS behavior of the homogeneous 
XY model with open boundary conditions (OBC), which 
is formally equivalent to the limit p — > oo of model p05p 
with L = 2l t + 1 and It integer. Note that translation 
invariance is lost in systems with OBC. Fig. [2] shows the 
data for x = (r, 0), which clearly support the expected 




FIG. 3: TSS plot of the two-point function G(0, x) at T c for 
a harmonic potential V = (r/lt) 2 . We set k = (top) and 
k = 1 (bottom). To check finite-size effects, we report two 
data sets: L m \l t and L ~ Sit in the two cases, respectively. 



FSS behavior 

G(0,x) = L- 1 / 4 (lnL) 1 / 8 .g(a;/i). (107) 

Note that OBC breaks translation invariance, and gives 
rise to power-law boundary corrections [ill |47| which are 
expected to be 0(L _1 ). In Figs. [3] and 0] we show the re- 
sults of the simulations for two different values of p, p = 2 
and p = 4. In order to check the scaling behavior (|106p . 
we present TSS plots with k = 0, using the analog of the 
FSS formula ()107|) . and with an optimal nonzero value of 
the coefficient k, which is determined by looking for the 
best collapse of the data. Optimal scaling is obtained be 
setting k « 1 for p = 2 and K ps 1/2 for p — 4, with 
an uncertainty which we estimate to be approximately 
10%. This simple scaling test clearly favors a nonzero 
p-dependent value for k. Taking also into account that 
k = for p — > oo, the above numerical results hint at the 
simple formula n — 2/p. These results should be univer- 
sal, hence they also apply to the BKT transitions of gen- 
eral systems of 2D interacting bosonic particles trapped 
by an external space-dependent potential, such as those 
which have been investigated experimentally 04IH HH ■ 
Finally, in Fig. [S] we compare the scaling functions 
associated with G(0,x) for p — 2, p = 4, and for 
the homogeneous system with OBC. We plot G r (x) = 
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FIG. 4: TSS plot of the two-point function G(0, x) at T c for 
a trap with potential 1/ = (r/l t ) 4 . We set k = (top) and 
k — 1/2 (bottom). To check finite-size effects, we report data 
for L « 4Z t and L « 2k- 
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FIG. 5: Plot of G r (x) = 4 /4 (ln«i)" 1/8_K/4 G(0,x) vs. r r = 
r(lnit) rc /Zt for p — 2, p = 4, and for the homogeneous system 
with OBC (it corresponds to p — s- oo), using the data for the 
largest available trap sizes. We use k = 1, 1/2, for p = 2, 4 
and the OBC case, respectively. 

?y 4 (ln Zi) _1 / 8_K / 4 G(0, x) versus the scaling variable r r = 
r(ln lt) K /h for the largest available trap size or lattice, us- 
ing k — 2 /p. For small r r all curves run very close up 
to r r of order one. This behavior is not surprising, since 
for r/l t — > boundary effects should become irrelevant, 
hence all scaling functions should behave in the same 



manner. 



VII. CONCLUSIONS AND SUMMARY OF THE 
MAIN RESULTS 

We have investigated the general features of the RG 
flow at the BKT transition, providing a definite char- 
acterization of the asymptotic critical behavior includ- 
ing the universal multiplicative or subleading logarith- 
mic corrections, in different critical regimes: (i) in the 
infinite-volume critical disordered phase; (ii) in the FSS 
limit, both for T > T c and T = T c \ (hi) in the TSS limit, 
again both above and at the critical temperature. 

For this purpose, we exploit the mapping between the 
standard XY or Coulomb gas models and the SG model, 
whose RG flow around the renormalizable fixed point de- 
scribes the BKT transition [3-4] . To determine the RG 
flow, we first derive a simple canonical universal expres- 
sion for the /3-functions. Then, we determine the asymp- 
totic solutions of the RG equations in the different sit- 
uations mentioned above. For the TSS case, numerical 
simulations confirm the RG predictions. 

We summarize our main results: 

(a) We prove that an appropriate analytical redefini- 
tion of the couplings of the SG model, cf. Eq. (ITTj) . allows 
us to write the /3-functions as 

/3 u (u,v) = -uv, (108) 
P v {u 1 v) = -u 2 [l + vf{v 2 )] 1 

to all orders of perturbation theory. The universal func- 
tion f(x) can be expanded as f(x) = bo + bix + b2X 2 + .... 
The zeroth-order term bo = 3/2 can be obtained from 
the two- loop results of Refs. [HI, HH, while the next coef- 
ficient b\ can only be determined by means of a four-loop 
calculation in the SG model. In practice, our results ex- 
tend the knowledge of the RG flow of the SG model to 
three loops. Under the same redefinition of the couplings, 
the anomalous dimension 7(tt, v) of the fundamental field 
becomes 7 = — 1/4 + v/8 — v 2 /16 + g2ou 2 + . . . , where 
1720 is a nonuniversal constant which does not enter the 
universal scaling behavior of the two-point function. 

(b) In the high-temperature critical regime and in 
the infinite- volume limit, the RG flow implies the BKT 
asymptotic behavior (setting r = T/T c — 1) [H-IH 

U{r)=Xe c '^[l + 0{^)] , (109) 

for the infinite-volume correlation length, where X is a 
nonuniversal constant. In the case of RG invariant quan- 
tities R, such as the ratio of two different definitions of 
correlation lengths and the zero-momentum four-point 
renormalized coupling [331 ] . we have 

R(t) =R*+ 0(1/ In 2 Coo). (HO) 

Corrections of order l/ln£oo are absent. The susceptibil- 
ity, defined as the space-integral of the two-point function 
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of the fundamental field, behaves as 



X 



- A P 1 ^ 



1 



o 



(111) 



where A x is a nonuniversal amplitude. The correction 
term behaving as 1/ In ^ is universal. 

(c) We have studied the FSS behavior. In the 
high-temperature phase, for any observable 0, the ra- 
tio 0(L)/0(L — > oo) approaches a universal function 
^(L/^oo) with 0(1/ In 2 L) corrections. The approach to 
the L — » oo limit is not uniform in T. At T = T c further 
logarithms appear. Indeed, setting // = ln(i/A) where A 
is an appropriate nonuniversal length scale, we have 



/' - 



i ln/x 



+ 0(//" 2 ), (112) 



where i?^ and are universal, depending only on the 
shape of the systems and their boundary conditions. An 
analogous formula is obtained for any RG invariant quan- 
tity, such as the helicity modulus and the Binder cumu- 
lants. Moreover, the asymptotic behavior of the suscep- 
tibility reads 

xL- 7/i = A x \u+-hxii) [l + OQi- 1 )], (113) 

where A x is a nonuniversal constant. Note two impor- 
tant differences between the behavior at T c and in the 
high-temperature critical regime (b). For RG invariant 
quantities corrections decay as 1/ In L in the first case; in 
the second one corrections of order 1/ In ^ are instead 
absent. Second, the behavior of the susceptibility is char- 
acterized by a leading logarithmic term at T c , while no 
such term is present in the high-temperature case, see 

Eq. (np. 

(d) Finally, we consider BKT transitions in 2D inter- 
acting bosonic particles which are trapped within a lim- 
ited region of space by an external space-dependent force, 
which is a common feature of experiments with diluted 
atomic vapors (5lj and cold atoms in optical lattices [52j . 
We investigate how the BKT critical behavior is affected 
by the presence of the external space-dependent trapping 
potential V(r) — {r/k) p coupled to the particle (energy) 
density. We consider observables derived from the two- 
point function G(x, y) of the fundamental field describing 
the critical modes, such as the one-particle correlation 
function {b^ x b y ) of the BH model (|74[) . which describes 
trapped bosonic atoms in an optical lattice. The analysis 
of the RG flow shows that TSS at the BKT transitions is 
characterized by a trap exponent 9 = 1, with additional 
multiplicative logarithms at T — T c . For example, at T c 
the two-point function scales as 

G(0,x) =Z t - 1/4 (lnZ t ) 1/8+K/4 5 P [r(ln/ t )7^, (114) 

where k is a new exponent which arises from the analysis 
of the RG flow of the external potential and which is ex- 
pected to depend on the power p characterizing the trap 



potential. Of course, in the limit p — > oo we must have 
k — > 0, since we must recover the known FSS behavior 
of a homogeneous systems. The scaling equation Q114p 
implies also that the correlation length £ t of the critical 
modes behaves asymptotically as 



&~Z t (lnZ t )- 



(115) 



These results are supported by Monte Carlo simulations 
of a 2D XY model with a space-dependent potential cou- 
pled to the energy density. They provide a clear evi- 
dence of the multiplicative logarithms in Eqs. (| 1 14[) and 
(|115p and are numerically consistent with the conjecture 
k = 2/p. 

Acknowledegment: We thank Gianpaolo Cicogna and 
Martin Hasenbusch for useful discussions and correspon- 
dence. 



Appendix A: Canonical form of the BKT /3-functions 

In the following we prove that the SG /3-fuctions can 
be simplified by a redefinition of the couplings, reduc- 
ing them to the canonical form given by Eqs. (|17[) and 
(p~8|) . For a general discussion of the mathematical prob- 
lem of the reduction of coupled differential equations to 
canonical form, see, e.g., Rcfs. [63, 64]. 

To all orders in the couplings a and 5, the ^-functions 
of the SG model have the generic form 

l3 a = -2aS+ ]T b a . nm a n S m , (Al) 



n+m>2 



1 2 i u nsr, 

+ 2^ °8,nmOt 



(A2) 



n+m>2 



In the SG model the sign of a is irrelevant, which implies 
the symmetry relations 



a (a,6) = -0 a (-a,6), 



(A3) 
(A4) 



As a consequence, b a , nm — if n is even and bs t nm = 
if n is odd. Moreover, for a = the theory is free and 6 
docs not flow. Hence 



fr(a = 0,5)=0, 



(A5) 



which implies bs -nm = if n = 0. 

We wish now to prove that, by an analytic redefinition 
of the couplings, 



n u n v m 



a = a a . 10 u+ ^ 

n+m>2 



S = cis.iqv + 2J as, nm u n v m , 

n-\-m > 2 



(A6) 
(A7) 
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we can rewrite the /3-functions of the SG model in the 
form 

/3 u (u,v) = -uv, (A8) 
= -u 2 (l + Y.b k v 2k+1 \. 

\ k=0 / 

To prove the general result, we shall work by perturba- 
tive induction. We assume that we have already proved 
the result to order n — 1, i.e. that we redefined couplings 
so that 



and 



p u = -uv + Hk-i[u,v) 



(A9) 



k — n 

m-2 



Pv = -u 2 1 + ]T hy2k+1 + u * G k-2(u, v), 

\ fc=0 / k=n 

where m = [n/2\, and H k (u,v) and G k (u,v) are homo- 
geneous polynomials satisfying 

H k (Xu,Xv) = X k H k (u,v), 

G k (\u,\v) = \ k G k {u,v). (A10) 

Moreover, they are even functions of u. We now consider 
the change of variables 



v! = U + uA n - 2 {u, w), 

v' = v + B n -x(u, v), 



(All) 



where A n ^ 2 (u, v) and B n -i(u,v) are homogeneous poly- 
nomials even in u and satisfying 

A n - 2 (Xu,Xv) = X n - 2 A n - 2 (u,v), 
B n ^i(Xu, Xv) = A n_1 S„_i(u, w). (A12) 

We wish to show that, by a proper choice of A n ^ 2 and 
_B„_i, we can cancel all terms of order n except, if n 
is odd, the term of order u 2 v n ~ 2 in f3 v . Hence we can 
obtain 



p a >(u',v') = -u'v' + u' J2 Hk-i(u',t/), (A13) 

k=n+l 

(m'-2 '' 
i + b ^ ,2k+1 

oo 

+ u' 2 G k - 2 (u',v'), 

k =n+l 

with different H and G and ml = [(n + 1 ) / 2 J . 

Requiring all terms (except the one mentioned above) 
to cancel, we obtain the equations 

uvA n _ 2 + uS n _! - ™ 9{u ^- 2) (A14) 
ov 



2u A n - 2 — uv — - u — „ h (Alo) 



du 

n-2 = Rn 



dv 



-u 2 G n - 2 = R 



where R n — for n even, and R n = b( n _ 3 y 2 u 2 v n 2 if n 
is odd. Eq. (|A14[) gives us the function B n -\. 

dA n -2 2 dA n - 2 , A1 „> 

B n _ 1 =uv— 5 r-u — ^ H n -i- (A16) 

ou ov 

Substituting it in Eq. (|A15|) . we obtain 
a 2 d 2 A n _ 2 , n 3 d 2 A n _ 2 A d 2 A 



du 2 



2u i v- 



dudv 



2 4" 
' U 



OA n r 



+u(u 2 + v 2 ) ' + 2u 2 v 
ou 



dA n r 



dv 



dv 2 

- 2u 2 A n _ 2 - R„ 



2rt dH n ^i 2 dH n _i 
u G n - 2 + uv — 1- u 



du 



dv 



Now, we expand 

n-2 

A n - 2 (u, v) = ^2 a k u k v n ~ 



2-k 



(A17) 



(A18) 



fc=0 



with a k = for k odd. Then, we must show that we 
can find ao, a 2 , so that the following equations are 
satisfied. For k even, with k > 4 and k < n — 1, we must 
satisfy 

E k = -k 2 a k + (4 - 3k + 2k 2 + 2n- 2kn)a k - 2 
-(2 + k 2 + in + n 2 - k(3 + 2n))a fe _ 4 - g k = 0. 

Here is the coefficient of order u k v n ~ k in the expansion 
of the r.h.s. of Eq. (|A17|) . Moreover, we should satisfy 

E n = (4 - n)a„_ 2 - 2a„_ 4 - = n even, 
= 3(3 - n)a„_ 3 - 6a„_ 5 - g„_i = rt odd, 
E 2 = 2(3 — n)ao — 4a2 — g 2 = rt even. 

i?2 cannot be satisfied for n odd — hence the necessity for 
the coefficients b k , which are defined by 

&(n-3)/2 = E 2 = 2(3 - n)ao - 4a 2 - g 2 . (A19) 

We now redefine a2fc for A; > 1 as 



a 2k 



k 

jfef J-- 1 - V 2 



(-l) fe ao 



II + '-^ 



(A20) 



With this redefinition, we obtain 



Ei = 6(4 - n)c 2 - 16c 4 - g 4 = (A21) 
B fc = -fc 2 c fe + (4 - 3fc + 2fc 2 + 2n - 2kn)c k - 2 

-(2 + k 2 + 3n + n 2 - fc(3 + 2n))c k - 4 - g k = 0, 

where 6 < k < n — 1 . For n even we should also consider 

E 2 = -Ac 2 - g 2 = 0, (A22) 
E n = p n ao + (4 - n)c„_ 2 - 2c„_ 4 - g n = 0, 
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where 



We define 



2(-l)«/2 



r r(n/2-l/2), (A23) 



(2-n)(n/2-2)! 

while for n odd we also have 

E n -i = 3(3 - n)c„_ 3 - 6c„_ 5 - <?„_i = 0. (A24) 

The parameter 6(„_3)/2 is defined by (n odd) 

6(n-3)/a = ^2 = -4c 2 - 5 2 . (A25) 

For n even it is evident that all equations can be solved. 
E n can be satisfied by fixing ao, while k < n — 2, 
can be satisfied by fixing Cfc. For n odd, one parameter, 
a , is no longer present: this explains why we are not 
able to satisfy all equations and we need to introduce 
the parameter 6( n _3)/2- In practice, we can satisfy £4 by 
fixing C4 as a function of C2, Eq by fixing cq and so on. 
This allows us to solve all equations except E n . However, 
there is still one free parameter, C2. Substituting the ex- 
pressions of c„_2 and c n _4 as a function of C2, we obtain 
E n = CX1C2 + ot%, with 



Oil 



(_l)(»-3)/2 2 (7-n)/2( n _ 2 )|, 

(n/2-3/2)! ■ 



(A26) 



Since ot\ 7^ 0, also E n can be satisfied, concluding the 
proof. 



Appendix B: Asymptotic expansions 

We wish now to discuss the computation of the asymp- 
totic behavior of integrals of the form 



dw 



h(w) 



Q + F(w) ' 



(Bl) 



where b > 0, a < 0, h(w) and F(w) are analytic functions 



and F(w) 



for 



0. The integral / diverges 



as Q — > if h(0) 7^ 0. The leading behavior can be 
computed by approximating F(w) m w 2 : 



dw 



h(0) nh{0) 



Q + w 2 



Q 



(B2) 



To compute the next nonanalytic term in the expansion, 
we consider 



[Q + F{w)f 



(B3) 



If x = wj y/Q, we consider the expansion in powers of 
y/Q at fixed x: 



h(xVQ) 



lQ + F(xVQ)} 2 



= E Q n/2 9n(x). (B4) 



n=— 4 



G(x,Q) = Q n/2 9n{x), 



(B5) 



i.e. the sum of the terms that diverge as Q — > 0, and 

g{w) =KmG{w/y/Q,Q) (B6) 

where the limit is taken at fixed w. It is easy to convince 
oneself that g(w) gives the principal part of the Laurent 
series of h(w)/ F(w) 2 , so that h(w)/F(w) 2 —g(w) is finite 
for w — > 0. We can thus rewrite 



J = 



h{w) 



[Q + F(w)Y 



dw 



dwG(w/^/Q,Q) 



G(w/y/Q,Q) 



(B7) 



The first integral is finite as Q — > 0, hence it does not 
contribute to the singular part of J. We can thus limit 
ourselves to considering the second term which can be 
rewritten as 



J 



b/VQ~ 



dx 



+ 



1 /V 2 (x)+.g_ 1 (x) 



- 3 / 2 g- 4 (x) + Q- 1 g_ 3 (x) + 



(B8) 



For x — > 00, we have g- n (x) ~ x~" (it is easy to prove 
it, using the fact that the expansion (|B4I) is indeed in 
powers of xy/Q). This implies that 



b/VQ~ 



a/VQ~ 



dx g- n (x) 



dxg~ n (x)+0{Q {n - 1)/2 ) (B9) 



for n = 2, 3, 4. For n — 1 we must be a little more careful. 
Assume that g—i(x) ~ ff- 1.00/2; for x — > 00. Then, 

rb/VQ 



cfeg_i(x) 

b/VQ 

dx 



g-i(x) 



g-i,oc>x 



1 



b/VQ 



dx 



a/VQ 



(B10) 

ff-l,ooZ 
l + Z 2 



The first integral now decays as 1/x 2 , hence we can ex- 
tend the integration limits to ±00 with corrections of 
order y/Q\ the second can be computed exactly and is 
finite for Q — >• 0. Hence the singular part of J is given 
by 



J 



+00 



d.v 



Q-^ 2 g^(x) + Q- x g- 3 {x) + 



+ Q- l/2 g-2{x)+g-i{x)- 9 -^ 

1 + x z 



(Bll) 



Using again the fact that the expansion (|B4|) is in powers 
of x^/Q, we observe the g2n(x) is even under x — > —x, 
while <?2n+i(a0 is odd. We obtain finally 



J = J-3/2Q 



-3/2 



J- lf 2Q- 1/2 + 0(l), (B12) 



15 



with 

J-Z/2 = 
J-l/2 

I f we now write F(w) 
J2 n h n w n , we obtain 



+00 



dxg- 4 (x), 



OO 
+ OO 



dxg- 2 {x). 



(B13) 
(B14) 



J2 n>3 F n w n , h(w) 



J. 



3/2 



nh 
2 ' 



(B15) 



J-1/2 = Yg (15F 3 2 /i - 12F 4 /io - 12F 3 /ii + 8/i 2 ). 

From the expansion of J we can easily derive the expan- 
sion of I: 

I = I-1/2Q- 1 ' 2 +Io + h/2Q 1/2 + 0(Q), (B16) 
with = 2J-3/2 and I 1/2 = -2J_ 1/2 . 



Appendix C: The canonical form of the trap 
/3-function 

We wish now to prove that the trap /3-function can 
be rewritten as in Eq. (|86j) . As in App. [A] we work by 
perturbative induction. We assume that we have proved 
the result at order n — 1, i.e. that the /3-function has the 
form 

w (u, v, w) = wT(u, v) + w 2 ^2 Hk-i{u,v,w), (CI) 



where T(0, 0) = 1 and Hk(u, v, w) are homogeneous poly- 
nomials of order k, i.e. satisfy 



Hk(Xu, Xv, Xw) — X k Hk(u, v, w). 



Then, we perform the change of variables 



z = w + wG n -i(u, v, w), 



(C2) 



(C3) 



where Gk(u,v,w) is homogeneous polynomial of degree 
k. If we now compute the /3-function associated with z, 



P z (u,v,z) = 



(C4) 



we obtain 



/3 Z = zT{u,v) + z 2 



H n - 2 (u,v,z) + 



dG n 



dz 



(C5) 



where we neglect terms of order n + 1 in the variables. 
Hence, if we define 



G„-i(u, v, w) 



dxH n -2(u,v,x), (C6) 



where the integral is performed at fixed u and v, we can- 
cel all unwanted terms, proving the result. 
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